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We construct supersymmetric Q-balls and boson stars in (d + 1) dimensions. These non-topological 
solitons are solutions of a scalar field model with global (7(1) symmetry and a scalar field potential that 
. appears in gauge-mediated supersymmetry (SUSY) breaking in the minimal supersymmetric extension of the 

Standard Model (MSSM). We are interested in both the asymptotically flat as well as in the asymptotically 
Anti-de Sitter (AdS) solutions. In particular, we show that for our choice of the potential gravitating, 
asymptotically flat boson stars exist in (2 + 1) dimensions. We observe that the behaviour of the mass and 
charge of the asymptotically flat solutions at the approach of the maximal frequency depends strongly on the 
number of spatial dimensions. For the asymptotically AdS solutions, the model on the conformal boundary 
can be interpreted as describing d-dimensional condensates of scalar glueballs. 
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^O. 1 Introduction 

A number of non-linear field theories possess solitonic-like solutions. These have broad applications in many 
. 1 branches of physics and constitute localized, globally regular structures with finite energy. Topological solitons 
\ P] possess a conserved topological charge that is connected to the existence of non-contractible loops in the 
1 theory. Non-topological solitons [21 [3] on the other hand appear in models with symmetries and possess a locally 
conserved Noether current and a globally conserved Noether charge. An example of such a non-topological 
soliton is the Q-ball [4] and its generalization in curved space-time, the boson star [51 [BJ [HI M DZ3 ■ These are 
solutions of models with self-interacting complex scalar fields and the conserved Noether charge Q is then related 
to the global phase invariance of the theory and is directly proportional to the frequency of the harmonic time- 
H ! dependence. Q can e.g. be interpreted as particle number [2]. As such, these solutions have been constructed 
in (3 + l)-dimensional models with non-renormalizable $ 6 -potential [111 1121 113[ . but also in supersymmetric 
extensions to the Standard Model (SM) [13]. In the latter case, several scalar fields interact via complicated 
potentials. It was shown that cubic interaction terms that result from Yukawa couplings in the superpotential 
and supersymmetry (SUSY) breaking terms lead to the existence of Q-balls with non-vanishing baryon or 
lepton number or electric charge. These supersymmetric Q-balls have been considered as possible candidates 
for baryonic dark matter [15] and their astrophysical implications have been discussed [TBJ- In [17], these 
objects have been constructed numerically using the exact form of a scalar potential that results from gauge- 
mediated SUSY breaking. However, this potential is non-differentiable at the SUSY breaking scale. In [TH] 
a differentiable approximation of this potential was suggested and the properties of the corresponding Q-balls 
have been investigated. Most models of a quantum theory of gravity need more than (3 + 1) dimensions and 
as such, it is surely of interest to investigate the properties of soliton solutions in higher dimensions. The first 
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study of Q-balls in higher dimensional space-time has been done in |19j . However, only a linearized version of 
the Lagrangian and equations of motion depending only in zeroth order on the ratio between the typical energy 
scale and the Planck mass has been used. In this case, an analytical solution can be given, however, the model 
does not capture the non-linear phenomena such as e.g. the behaviour of the mass and charge at the maximal 
frequency. A similar study has been done in [20] for (2 + 1) dimensions. Q-balls and boson star solutions of 
the full system of coupled non- linear equations in (4 + l)-dimcnsional asymptotically flat space-time have been 
investigated in [21] . Interestingly, it was found that the behaviour of the mass and charge at the approach 
of the maximal possible frequency is different for d — 4 as compared to d = 3. This was related to a scaling 
behaviour of the solutions at this critical approach and different dimensions of the spatial integrals. Spinning 
generalisations of these solutions can also be constructed [11] [T2] |13j [21] [22] [23] . These solutions possess a 
quantised angular momentum that is an integer multiple of the Noether charge. 

Topological and non-topological solitons in Anti-de Sitter (AdS) space-time have been investigated inten- 
sively recently. The interest in these objects is related to the AdS/CFT correspondence [2H US] which states 
that a gravity theory in a d-dimensional Anti-de Sitter (AdS) space-time is equivalent to a Conformal Field 
Theory (CFT) on the (d— l)-dimensional boundary of AdS. Interestingly, this is a weak-strong coupling duality 
that can be used to describe strongly coupled Quantum Field Theories with the help of weakly coupled gravity 
theories. This has been applied to a modeling of high temperature superconductivity with the help of classical 
black hole and soliton solutions in AdS 26, 2T][28]. The basic models use a scalar field coupled to a U(l) gauge 
field and the observation that close to the horizon of the black hole the effective mass of the scalar field can 
become negative with masses below the Breitenlohner-Freedman bound [29] such that the scalar field becomes 
unstable and possesses a non-vanishing value on and close to the horizon of the black hole. When computing 
the conductivities it turns out that the formation of a scalar field on a charged black hole corresponds to a 
phase transition from a conductor to a superconductor. However, insulator/superconductor phase transitions 
also play an important role in high temperature superconductivity and as such models including solitons have 
been suggested that describe this phenomenon [30] [3T] [32] . The AdS soliton is related to the black hole by a 
double Wick rotation with one of the coordinates compactified to a circle and has originally been suggested to 
describe a confining vacuum in the dual gauge theory [33] [34] since it possesses a mass gap. For solutions with 
Ricci-flat horizons there is a phase transition between the AdS black hole and the AdS soliton [35] which was 
interpreted as a confining/deconfining phase transition in the dual gauge theory. Note that this is different for 
black holes in global AdS where the black hole decays to global AdS space-time when lowering the temperature 

m- 

In the limit of vanishing gauge coupling the soliton solutions correspond to planar boson stars in AdS space- 
time. Since the scalar field is uncharged the interpretation in terms of insulators/superconductors is difficult 
in this case. However, since the AdS/CFT correspondence connects strongly coupled CFTs to weakly coupled 
gravity theories the prototype example of a strongly coupled field theory comes to mind - Quantum Chromo- 
dynamics (QCD). As such the planar boson stars in AdS have been interpreted as Bose-Einstein condensates of 
glueballs. Glueballs are color-neutral bound states of gluons predicted by QCD and the scalar glueball (which 
is also the lightest possible glueball) is predicted to have a mass of 1-2 GeV (see e.g. [37] for an overview 
on experimental results). Since these glueballs appear due to non- linear interactions and as such cannot be 
described by a perturbative approach it is very difficult to make predictions within the framework of Quantum 
Field Theory. However, holographic methods have been applied to make predictions about glueball spectra (see 
e.g. [38] and reference therein). 

Non-spinning boson stars in (d + l)-dimensional AdS space-time have been studied before using a massive 
scalar field without self-interaction [39] and in (3 + 1) dimensions with an exponential self-interaction potential 
[30]. Spinning solutions in (2 + 1) and (3 + 1) dimensions have been constructed in [53] and [JT], respectively. 

In this paper, we are interested in Q-balls and boson stars in both asymptotically flat as well as asymptotically 
AdS space-time with (d+ 1) dimensions. We use an exponential scalar field potential already employed in [40] , 
We will consider first the asymptotically flat case generalising some of the results obtained in 21] to higher 
dimensions and then also consider asymptotically AdS solutions. 

Our paper is organised as follows: In Section 2 we give the model, equations of motion and boundary 
conditions. In Section 3, we present our numerical results and conclude in Section 4. The Appendix 1 and 2 
contain results on the generalisation of an exact solution first found in [41] to (d + 1) dimensions and on the 
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existence of asymptotically fiat boson stars in (2 + f ) dimensions, respectively. 



2 The model 

In the following we will study non-spinning Q-balls and boson stars in a (d + l)-dimensional Anti-de Sitter 
(AdS) space time. The action S reads 

where R is the Ricci scalar, Gd+i denotes the (d + l)-dimensional Newton's constant, A is the negative cosmo- 
logical constant related to the AdS radius £ by A = — d(d — l)/(2£ 2 ). The second term on the right hand side 
of (JTJ is the Gibbons-Hawking surface term [32] with h the induced metric and K the trace of the extrinsic 
curvature on the AdS boundary. C m is the matter Lagrangian given by 

C m = -d M ®d M $* - U{\*\) , M = 0,l,....,d, (2) 

where $ denotes a complex scalar field and we choose the metric to have mainly positive signature. J7(|$|) is 
the potential 

Um) = m 2 vL y fl - exp (-^)) ■ (3) 

This potential is motivated by supersymmetric extensions of the Standard Model [TTJ [TS]. Here ?7 SUS y is a 

parameter such that rfeusy ^ corresponds to the energy scale below which supersymmetry is broken, while m 
denotes the scalar boson mass. 

The coupled system of ordinary differential equations is then given by the Einstein equations 

Gmn + A9mn = &KG d+1 T MN , M,N = 0,1,.. ,d (4) 

with the energy-momentum tensor 

dC 



Tmn — 9mn£ — 2 



Q g MN 

\g KL {d K <i>*d L <P + d L $*d K <S>) + U($) 



dM$*d N § + d N §*d M § (5) 



□ - 4^ ) $ = • (6) 



— 9M N 

and the Klein-Gordon equation 

The matter Lagrangian C m @ is invariant under the global U(l) transformation 

$ -> $e lx . (7) 
As such the locally conserved Noether current j M , M — 0, 1, .., d associated to this symmetry is given by 

j M = --($* d M ®-$d M ®*) with j$ = 0. (8) 
The globally conserved Noether charge Q of the system then reads 

Q = -J#xV=gf. (9) 
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2.1 Ansatz and Equations 

For the metric we use the following Ansatz in spherical Schwarzschild-like coordinates 

ds 2 = -A 2 {r)N{r)dt 2 + j^dr 2 + r 2 dtf d _ x , (10) 

where 

2n(r) 2A 9 

and dfl 2 ,_ 1 is the line element of a (d — l)-dimensional unit sphere. Note that the gravitational constant is 
chosen such that the d = 2 case has no Newtonian limit [43]. As such the metric function n(r) is well-behaved 
also in the d = 2 case. Note that for weak and static gravitational fields gu ~ — (1 + 2i/'(r)), where ip(r) is the 
Newtonian potential. Then the behaviour of V'( r ) i n d — 2 which is ip(r) ~ — ln(r) would imply the divergence 
of n(r) if that limit would exist within our Ansatz. 

For the complex scalar field, we use a stationary Ansatz that contains a periodic dependence of the time- 
coordinate t: 

d>(t,r) =e^(r) , (12) 

where w is a constant and denotes the frequency. 

In order to be able to use dimensionless quantities we introduce the following rescalings 

r — > — , u — > muj , l—t&lm , (f> — > ?7susy0 , n^nlm d ~ 2 (13) 
m 

and find that the equations depend only on the dimensionless coupling constants 

k = 8TrG d+1 T] 2 usy = 8tt — ^ , (14) 

where Mpi^+i is the (d + l)-dimensional Planck mass. Note that with these rescalings the scalar boson mass 
me = m becomes equal to unity. In these rescaled variables and coupling constants the coupled system of 
non-linear ordinary differential equations reads 

»' = "^(^" + t^) + ^) , (15) 

A' = kt (Ap + ^i) , (16) 
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(17) 



These equations have to be solved numerically subject to appropriate boundary conditions. We want to construct 
globally regular solutions with finite energy. At the origin we hence require 

0'(O) = , n(0) = , (18) 

while we choose A(oo) = 1 (any other choice would just result in a rescaling of the time coordinate). Moreover, 
while the scalar field function falls of exponentially for A = with 



4=rexp(-y/l-u 2 r) + ... (19) 

it falls of power-law for A < with 



(r » 1) ~ 
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When solving the equations numerically, we will choose as fourth boundary condition 0(oo) = for A = 0, while 
for A < we will choose the fall-off given in (|20p . 0a is then a constant that has to be determined numerically 
and which can be interpreted via the AdS / CFT correspondence as the value of the condensate of glueballs in 
the dual theory living on the d-dimensional boundary of global AdS. 
The explicit expression for the Noether charge reads 



oo 

Q = ^/drr^. (21) 
^ T(d/2) J AN v ' 



M = / dr r d - x N(f>' 2 + — t- + U(<p) . (23) 



o 

For k ^ 0, we can determine the mass M from the behaviour of the metric function n(r) at infinity. This reads 

n(r> 1) =M + ni r 2A+d + .... , (22) 

where ni is a constant that depends on t. 

For k = we have A = 1 and n = 0. Then the mass M corresponds to the integral of the energy density Tq 
and reads 

2 ^ /2 7, d i (mn " 2 <t> 2 
— — / dr r d ' x NA' 2 -\ - 

r(d/2)7 V n 

o 

Note that this expression is perfectly finite in AdS. While N contains a term cx r 2 the fall-off of the scalar 
function <j> guarantees that M has a finite value. Using the expression for the charge Q (f2"Tj) . we can give a 
relation between the charge and the mass 

M = ujQ + J dr ^ (iV0 ' 2 + U W ■ (24) 

o 

For A = it was found |18j that in the "thin-wall" approximation which corresponds to lo ~ the mass and 
charge are related as follows 

^ ~ *±± (25) 
This has been used to draw conclusions on the stability of the Q-balls. 



3 Numerical results 

The solutions to the coupled system of nonlinear differential equations are only known numerically. We have 
solved these equations using the ODE solver COLSYS [44]. The solutions have relative errors on the order of 

icr 6 -icr 10 . 



3.1 Q-balls 

We have first studied the case K = 0. This corresponds to Q-balls in a Minkowski (A = 0) or AdS background 
(A < 0), respectively. In this case, the Einstein equations decouple from the system and n(r) = 0, while 
A(r) = 1. 

It is known that Q-balls in d — 3 exist on a limited interval of the frequency u) G [w m m : w max ] . For A = 
the mass and charge diverge at both boundaries [121 [13], while for A ^ this is still true at w m ; n , but now the 
two quantities tend to zero at oj ma x 40 . This is related to the fact that in the limit uj — > uj max the scalar field 
function tends to zero everywhere. In Minkowski space-time this still leads to an infinite value of the integral 
since the space-time is infinite, however AdS space-time acts as a confining box and as such the integral becomes 
zero. 
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3.1.1 A = 



Q-balls with an exponential interaction potential of the form © have been studied in [TBI HO] in (3 + 1) 
dimensions. Here, we extend these results to d ^ 3 and show that the behaviour of the mass M and the charge 
Q depend crucially on the space dimension d. 

Our results for the mass M and charge Q in dependence on u are given in Figfl]for d = 2,3,4,5,6. We 
observe that the mass and charge diverge at u = w max for d > 3, while for d = 2 they tend to finite values. 
Moreover, since the potential and effective potential U c s ■= u! 2 <fi 2 — U(4>) do not depend on d the arguments 
employed in can also be used here, such that u; max = 1 and u> m i n = do not depend on d. This is clearly 
see in FigJTJ 




Figure 1: The value of the mass M (left) and the charge Q (right) of the Q-balls in dependence on the frequency 
uj in Minkowski space-time (A = 0) and AdS space-time (A = —0.1) for different values of d. 

In order to get an idea about the stability of these objects we can compare the mass M with that of Q free 
scalar bosons of mass tob = m. Due to our rescalings m = 1 and the mass of Q free scalar boson is just equal to 
Q. Any solution with M < Q would hence be stable to decay into Q free bosons. Our results for d = 2,3, 4, 5, 6 
are shown in Figj2] In d = 3 it was found [12l[T3j|40] that there exist two solutions with different M for a given 
charge Q, one of which is stable to the decay into Q free bosons, while the other is unstable. Here we observe 
that this is also true for d > 3: a stable branch exists up to a maximal value of the charge and then extends 
backwards to form a second branch that for some critical value of the charge becomes unstable. We observe 
that the bigger d the bigger is the value of the charge at which the solutions become unstable. On the other 
hand, for d — 2 we find that the solutions are always stable with respect to the decay into Q free bosons. The 
fact that Q-balls are stable for small values of u was pointed out already in [18]. In this so-called "thin- wall 
limit" the Q-balls fulfill the relation M ~ ((d+ l)/d)u>Q. Since lu is small M < Q in this limit and the Q-balls 
are stable with respect to a decay to Q free bosons. On the other hand for the "thick- wall limit" it is more 
difficult to make statement, but the results in [T8] again agree with our numerical findings. 

3.1.2 A^O 

This case corresponds to Q-balls in a fixed AdS background and has been studied for the exponential potential 
and d = 3 in [40] . Our results are shown in Fig 12] for A = —0.1 and d = 2,3,4,5,6. Similar to d = 3 the 
value of the mass M and charge Q tend to infinity at u> m - m — independent of d. On the other hand, the 
mass and charge tend to zero at w max . Moreover, we find that w max decreases with increasing d having the 
largest value for d = 2. For A ^ and n — it is known that in the d = 3 case and a particular choice of 
potential exact solutions to the scalar field equation exist [41]. We show in Appendix 1 that this generalizes 
to d dimensions. While the potential necessary to obtain this result is not of the form chosen in this paper, 
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Figure 2: We show the value of the mass M of the Q-balls in dependence on their charge Q for different values 
of d in Minkowski space-time. The small subplots show the behaviour close to the minimal value of Q. 



however, our numerical results for w max agree quite well with the analytic expression given by w max = A/£. 
This is related to the fact that for to — > u max the function <p(r) tends to zero everywhere. Hence all higher order 
terms in the potential become negligible and the Ansatz made in |41) gives a good result. 

Our numerical results for cj max in dependence on A and d are shown in Fig 13] together with the value of A/£. 
As is apparent from this figure the analytical result agrees quite well with our numerical values. Moreover, 
we observe as expected from the analytical result that w max increases with decreasing A and decreases with 
increasing d. 

In Figf4]we show the value of the mass M in dependence on Q for d = 2, 3, 4, 5, 6 and A = —0.1. Very similar 
to d = 3 the Q-balls have mass M larger than Q and are hence unstable, while for sufficiently large Q they 
become stable with respect to this decay. The value of Q = Q CI n at which this transition happens depends on 
A and it was found that Q C rit increases with decreasing A [40] . We find that the number of spatial dimensions 
d also has an influence on the value of Q C rit- We find that Q cl -it increases with increasing d. 

It is also known that radially excited Q-ball solutions exist which possess a number k E N of zeros in the 
scalar field function. In Fig 5(a) we show the mass M as function of uj and as function of Q, respectively for 
A = —0.1 and d = 3, 4 and k = 0, 1, 2. We observe that for fixed d the value of w max increases with the increase 
of k. Moreover, the bigger k the bigger is the difference between u; max for d = 3 and d — 4. The dependence of 
the mass M on Q shown in Fig 5(b) indicates that all solutions with nodes are unstable to decay into Q free 
bosons. This is not surprising since these can be seen as excited Q-balls in AdS space-time. 

As pointed out in [31] the field theory on the boundary of AdS describes condensates of scalar glueballs. 
This was further investigated in [40], where Q-balls in (3 + l)-dimensional asymptotically global AdS have been 
studied. In Figj6]we show our results for different values of A and d. Apparently, the expectation value of the 
dual operator < O > 1/ ^ A , which corresponds to the value of the condensate of scalar glueballs decreases for 
increasing d when fixing 0(0). This is related to the fact that the scalar field can spread into more dimensions 
when d is increased and hence less condensate is collected. Furthermore, the value of the condensate increases 
with decreasing A. This is connected to the fact that the "AdS box" decreases in size for decreasing A and as 
such the value of the condensate becomes bigger. 
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(a) w m ax over A 
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(b) Umax over d 



Figure 3: The value of w max in dependence on A (left) and in dependence on d (right). Though we plot d here 
as a continuous parameter, we should only read of the value for d £ N. We also give the value of A/£ and find 
that it gives a good approximation to our numerical data. 




Figure 4: We show the value of the mass M of the Q-balls in dependence on their charge Q for different values 
of d in AdS space-time with A = —0.1. 



3.2 Boson stars 

We now discuss the case k ^ 0. This corresponds to boson stars in an asymptotically flat (A = 0) or asymptot- 
ically AdS space-time (A < 0), respectively. 

3.2.1 A = 

As pointed out in [39], boson stars in (2 + l)-dimensional, asymptotically flat space-time do not exist for 
massive scalar fields without self-interaction. In the Appendix 2 we show that this is different in our case and 
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co Q 

(a) M over u) (b) M over Q 



Figure 5: The value of the mass M of the Q-balls in dependence on u> (left) and in dependence on the charge 
Q (right) in AdS space-time for different values of d and number of nodes k of the scalar field function. 




Figure 6: We show the expectation value of the dual operator on the AdS boundary < O > : / A corresponding 
to the value of the condensate of scalar glueballs in dependence on </>(0) for different values of A and d. 

that gravitating, asymptotically fiat boson star solutions in (2 + 1) do exist in our model with an exponential 
self-interaction potential. 

We have also studied the dependence of the mass and charge on the frequency u>. Our results are show in 
FigJ7]for the mass. The curves look qualitatively similar for the charge Q, this is why we don't show them here. 
We observe that the behaviour at w max depends crucially on the number of spatial dimensions d. For d = 3 
the mass and charge tend to zero, while for d — 4 they tend to a finite value. This has already been observed 
before and is confirmed with our type of potential. For d = 5 we find that now the mass and charge tend to 
infinity at the approach of w max . We have integrated up to values of the mass and charge of 10 7 and believe that 
approaching w max even closer, these values would further increase. This can be understood using the argument 
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Figure 7: The value of the mass M of the boson stars in dependence on the frequency to for A = and different 
values of d and k. The small subfigures show the behaviour of M at the approach of w max for d — 3, 4, 5 (from 
top to bottom). Note that the curves for the charge Q look qualitatively very similar, this is why we don't give 
them here. 

employed for d = 4 in |21] , As noticed in this latter paper, the scalar field function and radial coordinate show 
a scaling behaviour that is equal in d = 3 and d = 4. We find that this is also the case here and generalizes to 
d > 4 such that the behaviour is <f)(r) —> (j)(f), n —> 0, A —> 1 for uj — > w max with 



with 4>q a constant that tends to zero in the limit lj — > w max . Now this implies e.g. for the charge Q (the 
argument works similarly for the mass M): 



For d < 3 this tends obviously to zero, for d = 4 this becomes constant, and for d > 5 this tends to infinity, 
respectively for (f>o — > 0. In addition to this we observe that the approach to w max is not smooth in d — 5. This 
is shown in Fig. |8] For d = 3 and d = 4 the mass tends smoothly to zero and a finite value, respectively. For 
d = 5 we observe that the mass tends to a finite value on a lower branch of solutions, but that close to w max 
new branches of solutions exist. These are quite small in extend and in fact are barely noticeable for d > 5. As 
such, a second branch of solutions extends backwards from w max down to a critical value of co and then bends 
backwards to tend to infinity. Our conclusion hence is that while these solutions can exist for arbitrarily large 
values of the mass (and charge) there exists a mass gap in which solutions are not allowed. Furthermore, there 
is a small interval of u> in which up to three solutions with different masses exist. To understand this pattern, we 
plot the three solutions for lj close to w ma x in FigJS] We observe that for the same value of ui the three solutions 
are distinguished by the value of </>(0) with (f)(0) decreasing from the first to the third branch. Moreover, the 
solution spreads out more and more over r. On the first branch, the solution is still quite localized around the 
origin, while it becomes very delocalised on the third branch. 




(26) 




(27) 



o 



o 
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Figure 8: We show the value of the mass M of the boson stars in asymptotically flat space-time (A = 0) in 
dependence on the frequency ui close to iv max . Note that only d = 4 and d — 5 are physical values, but that 
the dimension is a parameter in our numerical programme that can also have non-integer values. Here we 
demonstrate how the mass evolves when going from d = 4 to d = 5. The subplot further demonstrates that 
there exists a mass gap in d = 5. 




200 400 600 800 1000 

r 



Figure 9: We show the profiles of the scalar field function </>(r)/</>(0) for the case where three branches of 
solutions exist close to oj max in d = 5. Here k = 0.001. 



We further observe that for all d the solutions exist down to a minimal value of the frequency w m i n - At 
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Figure 10: We show the value of the mass M of the boson stars in asymptotically fiat space-time (A = 0) in 
dependence on their charge Q for different values of n and d. 




Figure 11: We show the value of the mass Mj (uQ) of the boson stars in asymptotically flat space-time (A = 0) 
in dependence on the frequency u for different values of k and d. 



this value of uj, a second branch of solutions exists that has lower mass and charge for a fixed value of uj. This 
second branch extends back to a critical value of uj and then forms another branch. This spiraling has been 
observed before in d = 3 and is apparently also present in d ^ 3. 

We can also read off the dependence of uj m - m and uj max on d and k from FigJT] As is apparent, w max = 1 does 
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neither depend on k nor on d. On the other hand, the value of the minimal frequency w m ; n depends strongly on 
k and d. For fixed d it increases with increasing k, i.e. the stronger the interaction between the gravitational 
field and the scalar field the larger we have to choose the value of the frequency to find solutions. This is true 
for all d that we have studied. We also notice that the decrease in uj is bigger in higher d when increasing 
k by the same amount. For fixed k our results indicate that w m i n increases with increasing d. Normally, we 
would expect gravity in higher d to become weaker since it can leak into extra dimensions, however, here we 
are keeping Gd+i, i.e. the strength of the gravitational interaction constant. 

We show the dependence of M on Q for different values of k and d in FigfTU] As expected boson stars 
in asymptotically flat space-time are stable to decay into Q free bosons since these objects are gravitationally 
bound. This is also true for the case d > 5, where no bound on the mass M and charge Q exist. As is clearly 
seen from FigfTU] we find that also for d = 5 the curve is always below M = Q. 

For k = it was found that for small values of u> there is a relation between M and Q that depends on to 
and d [IS]. In FigfTTJ we plot M/(u)Q) as function of u> and find that only for d = 2 the approximation of flat 
space-time M/(toQ) ~ (d + l)/d is a good approximation for small U). For larger values of d the solutions do 
not exist for small U) and M/{ujQ) is always smaller than one. 

3.2.2 A^O 

We have also studied boson stars in AdS space-time. Our results for the mass M and charge Q in dependence 
on the frequency u are shown in Fig |12l for A = —0.1. 

o 

O I : 1 I : 1 



O - 

o 




(a) M over ui (b) Q over ui 



Figure 12: The value of the mass M (left) and the charge Q (right) of the boson stars in dependence on the 
frequency u) in asymptotically AdS space-time (A = —0.1) for different values of d and n. 

As for k = 0, we find that the mass tends to zero at uj — > w max , while at w -> w m in we find the spiraling 
behaviour typical for boson star solutions. This is true for all d > 3. Hence, while boson stars in asymptotically 
flat space-time with d > 5 dimensions can have arbitrarily large masses and charges, their mass and charge are 
bounded from above in asymptotically AdS space-time. 

The value of u> max depends on A and d, but not on k. It increases with decreasing A and decreasing d. Again, 
this can be explained with the fact that the maximal value of w can be very well approximated by w max = A/£. 

In terms of the holographic interpretation our solutions describe scalar glueball condensates including back- 
reaction. Our results for the expectation value of the dual operator on the AdS boundary are given in Fig ll3l 

We observe that the value of the condensate at a fixed value of the mass M increases with increasing k, i.e. 
backreaction of the space-time. This is true for all d. Moreover, there exists a maximal possible value of the 
condensate at the maximal value of the mass corresponding to intermediate values of 0(0). 
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(a) < O > 1 '' A over 0(0) (b) < O > X / A over M 



Figure 13: We show the expectation value of the dual operator on the AdS boundary < O > X / A corresponding 
to the value of the condensate of scalar glueballs in dependence on ^>(0) (left) and in dependence on M (right) 
for different values of k and d with A = —0.1. 

4 Conclusions and Outlook 

In this paper, we have studied (d + l)-dimensional Q-balls and boson stars in (asymptotically) flat as well as 
AdS space-time. We have used an exponential self-interaction potential for the scalar field motivated from 
MSSMs. We find that the behaviour of the solutions in asymptotically flat space-time at the approach of the 
maximal frequency depends strongly on the number of spatial dimensions d such that for d = 2, 3, 4 boson stars 
exist only up to a maximal value of the mass and charge, while for d > 5 the mass and charge of these objects 
can become arbitrarily large. However, our numerical results also indicate that for d — 5 a mass gap for the 
solutions exists. On the other hand, the value of the maximal and minimal frequency does not depend on d 
in this case. We show that the existence of boson stars in (2 + l)-dimensional asymptotically flat space-time 
depends crucially on the choice of the scalar potential. For our choice of potential, (2 + l)-dimensional boson 
stars always exist. 

We also generalize the exact Q-ball solution in AdS space-time found in [41] to d dimensions and find that the 
conclusion about the maximal value of u) drawn from this solution gives a good approximation to our numerical 
results. As such the maximal value of the frequency U) decreases with decreasing AdS radius and with d. 

We can also draw conclusions about the stability of these objects with respect to the decay into Q free 
bosons. While there are stable as well as unstable Q-balls in Minkowski space-time for d > 3, they are always 
stable in d = 2. In AdS space-time Q-balls are stable for small values of the charge and unstable for large values 
of the charge. Boson stars in asymptotically flat space-time are always stable, while for asymptotically AdS we 
find that only for d — 2 and above a critical ui the solutions become stable. For all other d the solutions are 
stable. 

With view to the holographic interpretation spinning solutions in asymptotically AdS space-time play an 
important role. These solutions has been constructed recently in d = 3 |41j . However, the Ansatz was chosen 
such that the boundary theory possesses no rotation and hence describes static glueball condensates. It is surely 
of interest to generalize this to find a holographic description of scalar glueballs possessing angular momentum. 

Acknowledgments We gratefully acknowledge support within the framework of the DFG Research Train- 
ing Group 1620 Models of gravity. 



14 



References 



N.S. Manton and P.M. Sutcliffe, Topological solitons, Cambridge University Press, 2004. 
R. Friedberg, T. D. Lee and A. Sirlin, Phys. Rev. D 13 (1976) 2739. 
T. D. Lee and Y. Pang, Phys. Rep. 221 (1992), 251. 
S. R. Coleman, Nucl. Phys. B 262 (1985), 263. 

D. J. Kaup, Phys. Rev. 172 (1968), 1331. 

E. Mielke and F. E. Schunck, Proc. 8th Marcel Grossmann Meeting, Jerusalem, Israel, 22-27 Jun 1997, 
World Scientific (1999), 1607. 

R. Friedberg, T. D. Lee and Y. Pang, Phys. Rev. D 35 (1987), 3658. 
P. Jetzer, Phys. Rept. 220 (1992), 163. 

F. E. Schunck and E. Mielke, Class. Quant. Grav. 20 (2003) R31. 
F. E. Schunck and E. Mielke, Phys. Lett. A 249 (1998), 389. 
M.S. Volkov and E. Wohnert, Phys. Rev. D 66 (2002), 085003. 

B. Kleihaus, J. Kunz and M. List, Phys. Rev. D 72 (2005), 064002. 

B. Kleihaus, J. Kunz, M. List and I. Schaffer, Phys. Rev. D 77 (2008), 064025. 

A. Kusenko, Phys. Lett. B 404 (1997), 285; Phys. Lett. B 405 (1997), 108. 



[1] 
[2] 
[3] 
W 
[5] 
[0] 

[7 

[8 

[9 

[10 

[11 
[12 

[13 

[14 
[15 

[16] K. Enqvist and J. McDonald, Phys. Lett. B 425 (1998), 309; S. Kasuya and M. Kawasaki, Phys. Rev. D 
61 (2000), 041301; A. Kusenko and P. J. Steinhardt, Phys. Rev. Lett. 87 (2001), 141301; T. Multamaki 
and I. Vilja, Phys. Lett. B 535 (2002), 170; M. Fujii and K. Hamaguchi, Phys. Lett. B 525 (2002), 143; 
M. Postma, Phys. Rev. D 65 (2002), 085035; K. Enqvist, et al, Phys. Lett. B 526 (2002), 9; M. Kawasaki, 
F. Takahashi and M. Yamaguchi, Phys. Rev. D 66 (2002), 043516; A. Kusenko, L. Loveridge and M. Sha- 
poshnikov, Phys. Rev. D 72 (2005), 025015; Y. Takenaga et al. [Super-Kamiokande Collaboration], Phys. 
Lett. B 647 (2007), 18; S. Kasuya and F. Takahashi, JCAP 11 (2007), 019. 

[17] L. Campanelli and M. Ruggieri, Phys. Rev. D 77 (2008), 043504; L. Campanelli and M. Ruggieri, 
Phys. Rev. D 80 (2009) 036006. 

[18] E. Copeland and M. Tsumagari, Phys.Rev. D 80 (2009) 025016. 

[19] A. Prikas, Phys. Rev. D 69 (2004) 125008. 

[20] A. Prikas, Nucl. Phys. B 690 (2004) 201. 

[21] B. Hartmann, B. Kleihaus, J. Kunz and M. List, Phys. Rev. D 82 (2010) 084022. 
[22] Y. Brihaye and B. Hartmann, Nonlinearity 21 (2008), 1937. 
[23] D. Astefanesei and E. Radu, Phys. Lett. B 587 (2004) 7. 

[24] see e.g. O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri and Y. Oz, Phys. Rept. 323 (2000) 
183; E. D'Hoker and D. Z. Freedman, Super symmetric gauge theories and the AdS/CFT correspondence, 
arXiv:hep-th/0201253 M. Benna and I. Klebanov, Gauge-string duality and some applications, arXiv: 
0803.1315 [hep-th]. 



see e.g. A. Kusenko, hep-ph/0009089 



15 



[25] J. Maldacena, Adv. Theo. Math. Phys. 2 (1998) 231; Int. J. Theor. Phys. 38 (1999) 1113. 
[26] S. S. Gubser, Phys. Rev. D 78 (2008) 065034. 

[27] S. A. Hartnoll, C. P. Herzog and G. T. Horowitz, Phys. Rev. Lett. 101 (2008) 031601; JHEP 0812 (2008) 
015; G. T. Horowitz and M. M. Roberts, Phys. Rev. D 78 (2008) 126008. 

[28] for reviews see C. P. Herzog, J. Phys. A 42 (2009) 343001; S. A. Hartnoll, Class. Quant. Grav. 26 (2009) 
224002; G. Horowitz, Introduction to holographic superconductors, arXiv:1002.1722 

[29] P. Breitenlohner and D. Z. Freedman, Annals Phys. 144 (1982) 249. 

[30] T. Nishioka, S. Ryu and T. Takayanagi, JHEP 03 (2010) 131. 

[31] G. T. Horowitz and B. Way, JHEP 1011 (2010) 011 jarXiv: 1007.37141 [hep-th]]. 

[32] Y. Brihaye and B. Hartmann, Phys. Rev. D 83, 126008 (2011). 

[33] E. Witten, Adv. Theor. Math. Phys. 2 (1998) 505. 

[34] G. Horowitz and R. Myers, Phys. Rev. D 59 (1998) 026005. 

[35] S. Surya, K. Schlcich and D. Witt, Phys. Rev. Lett. 86 (2001) 5231. 

[36] S. Hawkinga and D. Page, Commun. Math. Phys. 87 (1983) 577. 

[37] V. Crede and C. A. Meyer, Prog. Part. Nucl. Phys. 63 (2009) 74. 

[38] V. Mathieu, N. Kochelev and V. Vento, Int. J. Mod. Phys. E 18 (2009) 1. 

[39] D. Astefanesei and E. Radu, Nucl. Phys. B 665 (2003) 594 |gr-qc/030913l] . 

[40] B. Hartmann and J. Riedel, Glueball condensates as holographic duals of super symmetric Q-balls and boson 
stars, larXiv:1204.62"39l [hep-th]. 

[41] E. Radu and B. Subagyo, Spinning scalar solitons in anti-de Sitter spacetime, arXiv:1207.3715 [gr-qc]. 

[42] G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15 (1977) 2752. 

[43] N. J. Cornish and N. E. Frankel, Phys. Rev. D 43 (1990) 2555. 

[44] U. Ascher, J. Christiansen and R. D. Russell, Math. Comput. 33 (1979), 659; ACM Trans. Math. Softw. 7 
(1981), 209. 

5 Appendix 1: An exact solution and the value of oj max 

In analogy to the d — 3 case [H] we find that there exists an exact solution when choosing the scalar field 
potential of the form U(4>) = /J, 2 <j> 2 — X(j> 2k and n — such that 



where k = 1 + 1/A. In the limit w — > w m ax it is known that the function <j)(r) spreads out and becomes zero 
everywhere. Hence, we can read of the value of w max from (|28|) which gives 




(28) 



A/£ . 



(29) 
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6 Appendix 2: Existence of gravitating, asymptotically flat boson 
stars in 2 + 1 dimensions 

We follow the calculations in |39j and employ a scaling argument. For this, we use the reduced action which 
reads 



S= dr 



1 -n'A - r d ^A Ncj>' 2 - v 



\87rGd4 



NA 2 



■U(4) 



(30) 



where the prime denotes the derivative with respect to r. We now rescale the variable r — > fj,r and require that 
(dS/dA0U=i = 0. This gives 



dr r d - x A 



(d-2) 1- 



2Ar 2 



0' 2 + d • [/(</>) 



(d-l)(d-2) 
u 

The equation for <f) (jTTJ) can be rewritten as follows 



d- 



4(d- l)n 2(d- 2)Ar 2 \ w 



d(d- 1) / 



,2^2 



= (31) 



Integrating this with the help of the boundary conditions (|18I) .([T5 1) we find that 



(32) 



dr r d - L A N((j>y + 



ldU 



w 2 <j> 2 \ 



d<j> NA 2 ) 



= 



(33) 



Combining with (|31[) we get 

oo 

d ■ n 



dr r 



d-l 



1 



„d~2 



ll2 dfidu, TTr± , 

6 +2 2^-^ 



^^e 25 ( 2Ar 



iV 2 \d{d-l) 



(d-2)r 

T.d-2 



. 



(34) 



For A = and d — 2 this reads 



/ 



dr rA 



(35) 



Now A^t/)' 2 is always positive. In [39] it was argued that for U(<f>) = m 2 cj) 2 no gravitating asymptotically flat 
boson stars in d = 2 exist. This can be immediately seen when realizing that in this case the potential term 
vanishes and the integral can never become zero. However for our case with U(<fi) = 1 — cxp(— </> 2 ) the potential 
term gives 

/' I ;it~ 

(36) 



(^0-^))=exp(-0 2 )(0 2 + l) 



1 



This expression is negative definite (and vanishes for = 0) such that we can have gravitating asymptotically flat 
boson stars in d — 2. Note that for the 4> 6 self- interaction potential the negative-definiteness of this expression 
depends crucially on the choice of potential parameters. As such, our choice of potential also seems "more 
natural" from this point of view, since we can have solutions in all possible limits and dimensions of our model 
without "fine-tuning" the potential. 
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